Vortex quasi-crystals in mesoscopic superconducting samples 
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There seems to be a one to one correspondence between the phases of atomic and molecular 
matter (AMOM) and vortex matter (VM) in superfluids and superconductors. Crystals, liquids 
and glasses have been experimentally observed in both AMOM and VM. Here, we propose a vortex 
quasi-crystal state which can be stabilized due to boundary and surface energy effects for samples of 
special shapes and sizes. For finite sized pentagonal samples, it is proposed that a phase transition 
between a vortex crystal and a vortex quasi-crystal occurs as a function of magnetic field and 
temperature as the sample size is reduced. 

PACS numbers: 74.25.Op, 74.25.-q, 74.25.Dw 
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I. INTRODUCTION 

The general subject of vortex physics in superconduc- 
tors is quite interesting since there seems to be a large 
variety of possible equilibrium vortex phases in supercon- 
ductorsi. The term "vortex matter" has been coined to 
emphasize the complexity and diversity of vortex phases 
in superconductors when compared to atomic and molec- 
ular matter. One can think of a one to one correspon- 
dence between phases in atomic and molecular matter 
(AMOM) and phases in vortex matter (VM). A liquid 
in AMOM corresponds to a vortex liquid in VM2£; a 
crystalline lattice in AMOM corresponds to a vortex lat- 
tice in VM^; an amorphous or glassy solid in AMOM 
corresponds to an amorphous or glassy vortex system in 
VM£. In addition, a quasi-crystal is another very inter- 
esting state which has been experimentally discovered in 
AMOM£, but there are no corresponding experiments for 
vortex matter. 

The possibility of vortex quasi-crystals was initially 
discussed several years ago 7 -, but more recently quasi- 
periodic arrangements of vortices were discussed in 
the presence of underlying quasiperiodic pinning poten- 
tial^i-— — — . In these recent theoretical and experimental 
studies it was suggested that the physical properties are 
quite unconventional due to the quasi-periodic pinning 
potentials. However, unlike these studies, the present pa- 
per is dedicated to a discussion of vortex quasi-crystalline 
phases in superconductors without quasi-periodic pin- 
ning potentials. This manuscript is an extension of our 
previous unpublished workii£, where it was argued that 
vortex quasi-crystals maybe stabilized by boundary ef- 
fects and surface energies in finite systems. The present 
situation is quite distinct from the case of quasi-periodic 
pinning of vortices 8 - and from the case encountered in 
some atomic and molecular matter, where certain types 
of interactions can lead to colloidal quasi-crystals^. 

The central question of this manuscript is: under what 
conditions a quasi-crystalline arrangement of vortices is 
at all possible? A definite possibility is to argue that 
a stable vortex quasi-crystal can arise from an imposed 
quasi-periodic potential like for instance in the case of a 
superconductor with quasi-periodic pinning potentials 8 , 



or in heterostructures consisting of a superconductor film 
grown on top of a quasi-crystal film. Another possibil- 
ity is to create a quasi-periodic optical lattice with laser 
beams, trap and cool atoms, and produce vortex quasi- 
crystals in Bose or Fermi superfluids. This experiment is 
natural since the production of vortex lattices in super- 
fluid Bose or Fermi ultra-cold atoms has become stan- 
dar d and more recently square optical lattices were 
used to induce transitions between triangular and square 
vortex lattices using a mask technique^. Thus, the gen- 
eration of five-fold quasi-periodic optical lattices using a 
mask with five holes as the vertices of a regular pentagon 
(or other methods) may also be used to produce transi- 
tions between triangular and five-fold symmetric vortex 
lattices. 

However, in this manuscript, we concentrate on the 
possibility of stabilizing vortex quasi-crystals only due 
to boundary effects in finite superconducting samples, 
where the sample size and shape play an important role. 
This option is motivated by recent experimental studies 
of the vortex structure in disk, triangular, square, and 
star-shaped mesoscopic sample o 17 ' 18 i 19 ' 20 i 21 . These works 
provide us with the technology that allows the prepa- 
ration of pentagonal (pentagon-cylinders) or decagonal 
(decagonal-cylinders) samples as potential candidates to 
produce 5-fold vortex quasi-crystals. Furthermore, in 
these mesoscopic systems of cross-sectional area A, the 
number of vortices N is essentially given by N — HA/Qq, 
where $o is the quantum of flux, and H is the magnetic 
field. Therefore samples with large upper critical field 
H C2 may produce a large enough number of vortices even 
the cross-sectional area is small, and the appearance of 
quasi-crystalline order is a definite possibility. For exam- 
ple, for a sample of cross sectional area A = l/im 2 and 
upper critical field at zero temperature H C2 (0) = 10T, 
the number of vortices N 5000 is quite substantial close 
to H C2 (0). In addition, since superconductors with large 
H C2 are usually associated with short coherence lengths, 
the vortex quasi-crystalline phase proposed is more likely 
to be observable in proper samples of short coherence 
length superconductors at high magnetic fields and low 
temperatures. 

This possibility is considered here under the following 
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program. In this manuscript only the case of an isotropic 
type II superconductor in a magnetic field is considered. 
First, the bulk free energy is calculated for a triangular, a 
square and a 5-fold quasi-crystal array. The 5-fold quasi- 
crystal array is modeled by a Penrose tiling of the plane 
as shown in Fig. 1. It is shown that the Penrose tile 
array (vortex quasi-crystal) has a bulk free energy which 
is just a few percent higher than the triangular array. 
Second, instead of considering an infinite (bulk) system, 
a pentagon cylinder sample is discussed. The pentagon 
cylinder has a pentagonal cross-section (in the xy plane) 
with side dimension I, and with height L along the z di- 
rection. In this case, when the sample size gets smaller 
the contribution of the boundaries (surface energy) to the 
total free energy of the system becomes more important. 
The surface energy is highly sensitive to the symmetry 
and to the surface area of the boundaries. Taking into 
account the surface free energies, it is shown that the 
Penrose tiling (vortex quasi-crystal) has lower free en- 
ergy than the triangular lattice in certain regions of the 
magnetic field versus temperature phase diagram. This 
is suggestive that a "first order phase transition" may oc- 
cur between the triangular lattice and the Penrose tiling 
(vortex quasi-crystal) 2 ^.. 

The remainder of this manuscript is organized as fol- 
lows. In section UH we discuss a Ginzburg-Landau ap- 
proach to obtain the free-energy of the triangular, square 
and Penrose-tiling in the bulk. In section ITTTl we present 
the effects of pentagonal boundaries on the free energy 
of the triangular lattice and Penrose-tiling, and obtain 
the critical magnetic field where the transition from tri- 
angular vortex lattice to a vortex quasi-crystal occurs. 
In section HVl we discuss thermodynamic quantities like 
the change in magnetization and entropy, while in sec- 
tion [Vj we discuss neutron scattering, Bitter decoration 
and scanning tunneling microscopy as possible probes to 
identify the quasi-crystalline structure. Lastly, in sec- 
tion [VT] we state our conclusions and summarize. 



II. GINZBURG-LANDAU THEORY 

The starting point of our analysis to describe the 
vortex-quasicrystal state is the Ginzburg-Landau (GL) 
free energy density 
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for a bulk isotropic superconductor with no disorder, 
where AF S = F s — F n is the free energy density dif- 
ference between the superconducting state (F s ) and the 
normal state (F n ), and F x = a\^(r)\ 2 + /3|*(r)| 4 /2. 
It is useful to introduce the dimensionless quantities 
/ = *y^37^, p = r/A(T), A = 2^(T)A/$ , and 
U = 2tv£(T)\(T)H/$ , where A is the penetration depth, 
£ is the coherence length, and $o = hc/2e is the unit 
flux. Here, / = foexp(i<fi), and A = Aq + V^/k with 




FIG. 1: We consider a vortex lattice in a pentagon supercon- 
ducting sample having Penrose tiling with 5-fold rotational 
symmetry. The tile consists of two types of lozenges, one with 
internal angles 36° and 144° (thin lozenge) and the other with 
internal angles 72° and 108° (thick lozenge). 



k = A/£ is the GL parameter. Notice that when /o = 1 
the system is fully superconducting, and when /o = the 
system is normal, thus /q < 1 always. Considering the 
minimization of free energy density with respect to A and 
it is easy to arrive at equations for the dimensionless 
functions TL and /q. The microscopic field is 



U{x,y) = k 



1 



H-H c 
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with H parallel to the z direction, and the equation for 
/o can be written as 



V 2 (log<?) + 2 K 2 = 0, 



(3) 



where g — f§ is a positive definite function. Notice that 
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k for FL = H C2 and /o = (g = 0). The most 



general solution of Eq. ([3]) has the form 2 ^ 

g{x, y) = exp [~k 2 {x 2 + y 2 )/2] exp [7(1, y)} , (4) 

where ^y(x,y) satisfies Laplace's equation V 2 7(x,y) = 0. 
This means that 7(2;, y) is a harmonic function excluding 
the locations of vortices, and can be expressed as the real 
part of an analytic function of z — x + iy. This observa- 
tion has very important consequences for the microscopic 
field profile Tt(x, y) of Eq. ([2]), which depends strongly on 
the structure of g{x,y). 

The bulk Gibbs free energy density is given by^ 



G S {H,T) = G n {H,T) 
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where the parameter f3 = (g 2 ) /(g) 2 is a geometrical factor 
independent of k. The notation (■ ■ • ) indicates average 
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over volume. It is important to notice that (3 > 1 no 
matter what is the form of g(x, y) because of the Schwartz 
inequality. In addition, notice that the Gibbs free energy 
above is a minimum, whenever (3 reaches its minimum 
value. 

For the purpose of calculating the parameter (3 and the 
free energies corresponding to different vortex configura- 
tions, the analytical structure of g(x,y) in the complex 
plane is used to rewrite it as 



g(z,z) = exp(-K 2 zz/2)|P(z)|, 



(6) 



where P(z) — J^Y\iLii z ~ z i)- Here, each z% corresponds 
to a zero of g in the complex plane, and M is the number 
of zeros. The zeros z< indicate the location of vortices. 
From now on it is assumed that there is only one vortex 
with flux $p at each position Zi, i.e., each zero is non- 
degenerate. In this case, M corresponds to the number 
of vortices, and thus the total flux threading the sample 
is $ = M$o. The normalization coefficient Af just guar- 
antees that g(z,z) < 1. Depending on the locations of 
the zeros of g(z, z), it is possible to study several possibil- 
ities of periodic and quasi-periodic vortex arrangements. 
In this study only vortex crystals corresponding to trian- 
gular and square lattices and vortex quasi-crystals corre- 
sponding to the 5-fold Penrose tiling of the plane are con- 
sidered. Both the square lattices and triangular lattices 
can be generated via the tiling method, i.e., via the pe- 
riodic arrangements of identical square tiles or identical 
lozenges of internal angles (60° and 120°). The Penrose 
lattice, however, requires quasi-periodic arrangements of 
two types of tiles (lozenges), one with internal angles 36° 
and 144° and the other with internal angles 72° and 108°. 
Using the representation in Eq. ((BJ) , the values of (3 for 
the triangular, square and Penrose tiling are respectively 
03 = 1.16, [3 A = 1.18 and f3 5 = 1.22. This immedi- 
ately indicates that the triangular lattice has lower free 
energy than the square lattice which has lower free en- 
ergy than the 5-fold vortex quasi-crystal (Penrose tiling) , 
as expected. However, the fact that the free energy dif- 
ference between the triangular and 5-fold vortex quasi- 
crystal is only a few percent suggests that appropriate 
boundaries can favor 5-fold symmetry as the sample size 
gets smaller, as can appropriately imposed quasi-periodic 
potentials. Thus next, we discuss boundary effects and 
the magnetic field versus temperature phase diagram. 



III. BOUNDARY EFFECTS AND PHASE 
DIAGRAM 

In order to investigate how boundary effects can mod- 
ify the total free energy of the system, a sample in the 
shape of a pentagon cylinder of side i and height L is 
considered. Imposing that no currents flow through the 
sample boundaries leads to the condition 



at all five side faces. The unit vector n points along the 
normal direction of each facet of the pentagon cylinder. 
The boundary conditions can be translated in terms of 
the harmonic function 7(2, y) as 
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where n x and n y are the x and y components of the 
normal unit vector n at each one of the pentagon cylinder 
side faces. The solution for this boundary value problem 
can be obtained using a Schwarz-Christoffel conformal 
map of the pentagon into a semi-infinite plane 



^=C(w- Xl )- 2 / 5 (w 2 
aw 
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where the vertices of the pentagon located at are 
mapped into the points (x±, ±%2, ±£3) on the real axis 
of the w-plane. The full solution of this problem is com- 
plicated, and requires heavy use of numerical methods 2 ^. 
However, the free energy density can be estimated when 
the bulk solution in Eq. ([6]) is treated as a variational 
solution of the boundary value problem determined by 
Eqs. ©, ® and ©. 

The Gibbs free energy density difference AG = G3 — 
G5 between the triangular and the 5-fold quasi-periodic 
Penrose structure then becomes 



AG 
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where (3* — (3^(3^/(13^ — ^3), H c is the thermodynamic 
critical field, and e = (a 3 — a^){2R e + 03 + a 5 )/2R 2 , 
with 013 = 0.93ao, «5 = 0.90eto, and ao = \/&o/H. The 
effective length of the sample R e = i/V^ — t 2 , where r = 
2cos(-7r/5) is the golden mean. Here, R e corresponds to 
the radius of the circle that circumscribes the pentagonal 
sample. 

This expression for AG is valid only when H 3> $0/ Rl- 
The second term in AG takes into account the bound- 
ary mismatch energy, and indicates that as the size of the 
pentagon cylinder gets smaller it becomes more favorable 
to have a 5-fold quasi-crystal rather than a regular tri- 
angular lattice. Notice, however, that when R e — > 00 the 
triangular lattice has lower Gibbs free energy as it must, 
and no transition to a 5-fold quasi-crystal occurs. Thus, 
this possible transition may occur for finite sized samples 
only. From the condition that AG = we obtain 



Ho — H r 



1 - K 
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n- [-iV - 2vrA/$ ] V(x,y) = 



(7) 



where the transition to a quasi-crystal occurs. Here, n* — 
\/2k 2 - 1/k. 

The phase diagram for a superconductor with k = 20, 
H C2 (0) = 10 T and R e = 10 -6 m is shown 25 in Fig. [2] 
using H C2 (T) = H C2 (0) [l - (T/T c ) 2 ]. Notice that the 
number of vortices N close to Hq can be quite large 
at low temperatures. For the parameters above N = 
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T/T c 

FIG. 2: H — T phase diagram for a pentagonal cylinder cut 
out of a regular cylinder of radius for R e — 10 _6 m. The solid 
and dashed lines represent H C2 and Hq, respectively, while 
H C1 is not shown. The superconductor is assumed to have 
K = 20, H C2 (0) = 10T. 



HqA/^q « 4000, and the vortex quasi-crystal structure 
can be extracted. However, the phase diagram becomes 
less precise at high temperatures, since the low value of 
H sets a restriction on the number of vortices. Further- 
more, for fixed H C2 (0), one can find the critical value of 
R e , c below which the vortex quasi-crystal phase appears. 
Curves of i? e)C versus T or H can be obtained and are 
shown in Fig. [3] Notice that the critical R ec increases 
with increasing H (T) for a fixed T (H). Notice also 
that for a fixed sample size, and fixed magnetic field the 
vortex-quasi crystal phase always occurs at higher tem- 
peratures due to its higher entropy. 

In order to address the signature of the crystal- 
quasicrystal phase transition, we discuss next thermo- 
dynamic properties including magnetization and entropy 
of the triangular and Penrose lattices in the pentagonal 
geometry. 
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FIG. 3: Critical value of superconducting sample size R £lC as 
a function of (a) H for a given T < T c (0), and as a function of 
(b) T for a given H < H C2 (0). Parameters used in these plots 
are the same as in Fig. [2j and the regions labeled TL, VQC, 
and N correspond to triangular lattice, vortex quasi-crystal, 
and normal phases, respectively. 



IV. THERMODYNAMICS: CHANGES IN 
MAGNETIZATION AND ENTROPY 

The jump of the magnetization AM(= M3 — M5) as 
a function of temperature at the critical field Hq can be 
calculated from the Gibbs free energy leading to 



An 



(2k 2 - 1)13* 



AM, 



(12) 



where 7* = K*y/JJ*e and AM S = (a 3 - 05) (Re + a 3 + 
a 5 )/An 2 Rl(l - 7*). A plot of AM is illustrated in Fig. 
[4] for the same parameters of Fig. [2] Using these pa- 
rameters produces jump discontinuities AM « — 0.060G 
at T = 0, and AM w -0.028G at T = 0.8T C . How- 
ever, measurements of AM may require the preparation 



of an array of identical pentagonal cylinders to enhance 
the overall value. Notice that AM < indicates that the 
5-fold vortex quasi-crystal is denser than the triangular 
vortex lattice at Hq, being at best a few percent denser 
at T = 0. 

In addition to magnetization measurements, it is also 
interesting to perform calorimetric experiments. How- 
ever, specific heat measurements are very difficult be- 
cause they require large samples. Since sample size is im- 
portant for the present discussion it is not clear that such 
experiments can be successfully performed. Nevertheless, 
the thermodynamic relationship between the magneti- 
zation and entropy jumps is revealed in the Clapeyron 
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FIG. 4: The jump discontinuity AM = M3-M5 at the critical 
field Hq for various reduced temperatures T/T c , using the 
same parameters of Fig. [2] 



equation 

AS = S 3 - S 5 = -AMdHg/dT. (13) 

Since dH Q /dT < 0, and AM < implies that AS < 
0, the entropy S3 of the triangular vortex lattice is less 
than the entropy 5s of the 5-fold vortex quasi-crystal, 
indicating that latent heat L = TAS is required to cause 
this phase transition. 

Thermodynamic quantities can provide a good un- 
derstanding of properties averaged over the entire sam- 
ple, and can characterize the signature of the crystal- 
quasicrystal phase transition. However, as discussed 
above, a good measurement of these quantities may re- 
quire an array of identical samples, which introduces ex- 
perimental complexity and difficulty. As another pos- 
sibility, the use of local probes discussed next is much 
desired in order to reveal the change in structure from a 
triangular vortex crystal to a 5-fold vortex quasi-crystal. 
For instance, neutron scattering, Bitter decoration or 
scanning tunneling microscopy (STM) experiments may 
help elucidate the structure of the vortex arrangement in 
mesoscopic samples. 

V. LOCAL PROBES: NEUTRON DIFFRACTION 
AND SCANNING TUNNELING MICROSCOPY 

In neutron diffraction experiments periodic or quasi- 
periodic variations of TL{x 1 y) will result in Bragg peaks. 
The position of these peaks determine the characteristic 
length scale of the vortex structure and its symmetry. 
The neutron scattering amplitude in the Born approxi- 
mation is 

M r 

&(q) = ^2 J »nH{v) exp(*q • r)dr, (14) 

where fj, n — 1.91eh/M n c is the neutron magnetic mo- 
ment and the M n is the neutron mass. The scattering 



amplitude 6(q) is directly proportional to the Fourier 
transform H(q) [H(q x ,q y )] of the microscopic field H(r) 
[H(x, y)] of Eq. ([2|). The neutron scattering cross section 

v(qx,q y ) = A-K 2 \b{q x ,q y )\ 2 (15) 

has sharp peaks at (q x , q y ) = (0, 0) (central peak) and at 
(qx,q v ) = (±qxNm,±q y Nm) (first Bragg peaks), where 
qxNm = Qn cos(mn/N) and q y Nm = Qn sin(m7r/iV), 
with m = 0, 1, N — 1. For the triangular lattice N = 3 
the first Bragg peak occurs at \Q^\ — 2.31 x ir/da, where 
c?3 is the lattice spacing. For the 5-fold vortex quasi- 
crystal (Penrose Lattice, N = 5) the first Bragg peak 
occurs at |Qs| = 2.46 x ir/d^, where d$ is the side of a tile. 
Since the sample size is important for the observation of a 
5-fold quasi-crystal, neutron scattering experiments may 
be difficult to perform. 

However, Bitter decoration might be a useful tech- 
nique if magnetic nanoparticles could be used to deco- 
rate the magnetic field profile, and then be seen by a 
scanning tunneling microscope (STM) (magnetic or non- 
magnetic). Furthermore, it may be possible to use just 
an STM to scan over the pentagonal sample and probe 
the local density of states which is substantially different 
inside and outside of vortex cores, due to the presence of 
vortex cores states. In this case, it may be also useful 
to make a periodic pattern of pentagonal samples, to ob- 
tain an ensemble average. It should be possible as well to 
perfom STM scans at different fields and temperatures in 
the vicinity of Hq(T), which would reveal the real space 
locations of vortices. The pattern obtained could then 
be Fourier transformed (FT) to obtain a 6-fold pattern 
for the triangular vortex lattice and a 10-fold pattern for 
the 5-fold vortex quasi-crystal (Penrose lattice), which is 
shown in Fig. [5] 
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FIG. 5: First peaks of the square of the Fourier transformed 
(FT) pattern for the 5-fold vortex quasi-crystal (Penrose lat- 
tice). Notice that the pattern is 10- fold symmetric, unlike 
the case of the triangular vortex lattice, where the pattern is 
6- fold symmetric. 

Now that the phase diagram, thermodynamics, and 
the local signatures of a vortex quasi-crystal have been 
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discussed, it is important to say a word on the stability 
of such structures. A stability analysis in the free en- 
ergy can be performed by moving the vortices away from 
their equilibrium positions Zi to Zi + 8zi . The eigenval- 
ues associated with these displacements indicates that for 
H > Hq(T) the vortex quasi-crystal lattice is stable, and 
therefore potentially realizable in mesoscopic samples. 

VI. FINAL COMMENTS AND CONCLUSIONS 

Before concluding, we would like to make several com- 
ments in connection with the approximations used and 
the observability of the vortex quasi-crystal phase dis- 
cussed. 

First, it should be emphasized that our free energy 
analysis provides a preliminary understanding of the vor- 
tex quasi-crystal phase, but further detailed numerical 
work is necessary. For instance, the vortex quasi-crystal 
sitting on a Penrose lattice is only one possible state in 
a pentagon superconducting sample. Although we have 
shown that this state has lower free energy than a tri- 
angular or rectangular lattice, we cannot rule out other 
possibilities based on the present calculation. 

Second, since there is no full compatibility of the 5-fold 
Penrose lattice with the pentagon cylinder geometry, the 
appearance of disclinations and dislocations is possible. 
As a result there is an additional possibility of a solid 
(crystal or quasi-crystal) vortex structure in the center of 
the sample, which melts or gets disordered at the sample 
boundaries 2 ^. 

Third, the number of vortices in superconducting sam- 
ples is restricted by N max = H C2 (0)A/Qq, which means 
that for superconductors with high H C2 a large num- 
ber of vortices is possible. One can make a simple es- 



timate for -/Vmax by using the Ginzburg-Landau relation 
H C2 = $o/(2tt£ 2 ), leading to A max ~ [R e /t,{T = 0)] 2 . 
For conventional type II superconductors, the number 
A max is only about 10 1 ~ 10 2 , thus the thermodynamic 
limit and hence the definition of the quasi-crystal phase 
is questionable. However, for materials with short coher- 
ence length (and large upper critical fields), such as high- 
T c superconductors, where £(T = 0) is about 10 1 ~ 10 2 A, 
then A max can be as large as 10 4 ~ 10 5 , and the quasi- 
crystalline structure can be well defined and observed. 
Thus, the observation of a vortex quasi-crystal state is 
more likely to occur in short coherence length supercon- 
ductors at low temperatures and high magnetic fields. 

Finally, strong disorder in the sample can also destroy 
the quasi-crystal structure due to the pinning of vortices, 
however, clean mesoscopic superconducting materials al- 
ready exist and shell structures have been observed for 
Niobium samples of \im sizes 2 ^. Thus, we suspect that 
experimentally this should not be an issue as vortex lat- 
tices (triangular) are routinely observed in reasonably 
clean superconducting samples. Therefore, the choice of 
pentagonal mesoscopic samples of superconductors with 
sufficiently large H C2 should allow for the observation of 
the vortex quasi-crystal state. 

In summary, we have shown that vortex quasi-crystals 
may be experimentally observed in mesoscopic samples 
of type II superconductors with large upper critical fields 
(short coherence lengths). By taking into account bound- 
ary effects, sample shape and size, we proposed that a 
first order phase transition occurs between a vortex crys- 
tal and a vortex quasicrystal, as magnetic field and tem- 
perature are varied. 

We would like to thank Wai Kwok and Franco Nori for 
references, the Aspen Center for Physics for their hospi- 
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